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Abstract 

We present the timing characteristics of the flash ADC readout of the GlueX forward calorimeter, which depends on 
precise measurement of arrival time of pulses from FEU 84-3 photomultiplier tubes to suppress backgrounds. The tests 
presented were performed using two different 250 MHz prototype flash ADC devices, one with eight-bit and one with 
twelve-bit sampling depth. All measured time resolutions were better than 1 ns, independent of signal size, which is 
the design goal for the GlueX forward calorimeter. For pulses with an amplitude of 100 mV the timing resolution is 
0.57 ± 0.18 ns, while for 500 mV pulses it is 0.24 ± 0.08 ns. 
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' 1. Introduction 

The forward electromagnetic calorimeter (FCAL) of 
the GlueX experiment jTj at Jefferson Lab will use flash 
analog-to-digital converters (FADC) to detect signals from 
final-state photons, produced in photon-nucleon reactions. 
The calorimeter will consist of 2800 4 cm x 4 cm x 45 cm 
type F8 lead glass blocks, each coupled to an FEU 84-3 
photomultiplier tube (PMT) . Electrons and positrons pro- 
duced in the electromagnetic shower from an incident pho- 
ton generate Cherenkov radiation that will be converted 
to a current pulse by the PMT. The PMT pulses will 
be digitized by twelve-bit 250 MHz multi-channel FADC 
boards [2] . Continuous digitization by the FADC preserves 
the signal pulse shape, enabling simultaneous measure- 
ment of total charge, which is the integral of the pulse, 
and signal arrival time. The integrated charge is propor- 
tional to the energy that the incident particle deposited in 
the block and the pulse shape can be used to determine 
the arrival time of the particle. 

The pulse arrival time will be used to (1) isolate final- 
state photons from out-of-time hits in the detector that re- 
sult from the expected high level of electromagnetic back- 
ground and (2) determine the exact beam bunch that pro- 
duced the collision of interest. The electromagnetic back- 
ground arises from the photon beam, which will be deliv- 
ered at an average rate of 10^ 7/s incident on the target 
(at energies of approximately 9 GeV) with beam bunch 
spacing of 2 ns. Predicted background rates on the FCAL 
modules closest to the beam axis are of order 1 MHz; else- 
where in the detector the background rate is at the order 



of 10 kHz. Given the bunch spacing, a timing resolution 
of better than 1 ns for each FCAL block is desired to both 
suppress electromagnetic background and perform bunch 
identification. Note that a typical event will illuminate 
multiple blocks, thereby providing numerous independent 
measurements of the event time, each with better than 1 ns 
resolution, which will allow a determination of the beam 
bunch that produced the event of interest. 

In order to make optimal use of the data bandwidth 
available, it is critical that the timing measurement of the 
pulses be done in real-time as part of the digital signal 
processing of the FADC buffer for each event. We will 
present two such processing algorithms to determine the 
timing information from a PMT signal, both of which uti- 
lize samples from the leading edge of the PMT pulse. We 
demonstrate that both of these algorithms provide time 
resolution sufficient for the needs of the GlueX experiment. 

2. Hardware characteristics 

The readout of the GlueX forward calorimeter will con- 
sist of 2800 Russian-made FEU 84-3 PMTs, each pow- 
ered by a 24 V Cockroft-WaUon base g]. The PMTs will 
be connected to multiple 16-channel 250 MHz twelve-bit 
FADC boards, designed and produced by Jefferson Lab [5]. 
For tests described in this article, we utilized a single- 
channel eight-bit 250 MHz FADC developed at Indiana 
University in addition to a prototype eight-channel twelve- 
bit 250 MHz FADC developed by Jefferson Lab. We will 
briefly describe these two FADC devices and the pulse 
shape of the FEU 84-3 phototube below. 
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2.1. Single- and multi-channel FADC 

The eight-bit single-channel FADC prototype has been 
implemented as a PCI card so that it can be tested us- 
ing a standard personal computer. The FADC digitizer 
samples the input signal at 4 ns intervals, and the samples 
are continuously written into a dual-port ring-buffer. This 
enables the extraction of the samples without interrupting 
the digitizer. The data are read out over the PCI bus. 
The three processes, digitization, extraction and read out, 
are independent of each other and can happen simultane- 
ously, which allows the device to have zero dead-time. The 
analog input saturates for voltages above wl.lS V. 

The twelve-bit 250 MHz FADC is based on the MAX1215 
chip with a 1.45 V dynamic range. Like the eight-bit 
board, this device is also a zero dead-time device. It has 
the ability to read signals simultaneously from up to eight 
PMTs, with the digitization of all channels synchronized 
to the same clock. The board is based on the VME64x 
platform. More details on the operation and characteris- 
tics can be found in Ref. [2|. 

2.2. Characteristics of the FEU 84-3 PMT pulse 
Figure [l] shows a typical pulse from a FEU 84-3 PMT 

obtained by digitizing the signal every 0.2 ns with a digital 
oscilloscope. The line in Fig. [l] represents the fit obtained 
using a sum of a bifurcated Gaussian and a constant _B, 
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where A represents the amplitude of the signal, tp is the 
peak arrival time and cti_2 are the Gaussian widths that 
model the duration of the rising and falling edges, respec- 
tively. It is evident that the rising edge is well reproduced 
by a Gaussian line-shape, and this feature will be employed 
by one algorithm for the online extraction of the pulse ar- 
rival time. In addition, pulse characteristics are given by 
the parameters of the bifurcated Gaussian fit. We found 
that the characteristic leading-edge time (cri) is about 9 ns 
independent of the pulse amplitude. 

Since the pulse varies slowly around the peak, the peak 
sample time is a poor choice to characterize the arrival 
time of the pulse. We choose a characteristic signal time 
corresponding to the time at which the pulse has reached 
half of the signal maximum. We call this time to and 
discuss below two algorithms to measure to- 

3. Timing algorithms 

We explored two methods to obtain the characteristic 
time of a pulse that both rely on two FADC samples on 
the rising edge of the pulse. Other methods have been 
used, e.g., using the first algebraic moment [J; however, 
we have had obtained the most precise results using the 
methods outlined below. In the discussion that follows we 
assume that all FADC measurements have been adjusted 
to remove the portion of the signal due to the DC offset, 
i.e., all FADC samples are "pedestal subtracted." 
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Figure 1: A typical pulse from FEU 84-3 PMT obtained from a 
digital oscilloscope. (The pulse has been inverted for comparison 
with FADC data.) The solid line represents a fit to a bifurcated 
Gaussian with a constant offset. 



3.1. Gaussian transformation 

This method takes advantage of the fact that the rising 
edge of the pulse from the PMT is Gaussian in shape. 
Initially proposed by Teige et al. 5j it is based on the 
transformation that turns a Gaussian edge into a straight 
line: 



^' = W-ln 



(2) 



where Sp is the peak sample, while S'^ is the transformed 
sample obtained from an original sample 5*^ at the time ti . 
The transformed rising edge samples are now a linear 
function of time 



S'i = act J -I- be 



(3) 



with the slope, ac, and the intercept, bQ, related to the 
parameters of the rising-edge Gaussian, Eq. [T] The pa- 
rameters of the straight line can be calculated using any 
two samples from the rising edge. The line, in turn, deter- 
mines the time when the pulse reached any given fraction 
of its peak value. In particular, tg, is given by 



in = 



^fHf) 



ac 



(4) 



It was found that the two samples immediately preced- 
ing the peak sample give the best measure of to in this 
approach [S]. 

3.2. Linear interpolation 

It was desirable to test a different algorithm which 
did not involve mathematical operations such as a square- 
root and a logarithm in order to facilitate adaptation of 
the algorithm to the FADC field-programmable gate-array 



(FPGA). For that reason we considered another method 
that assumes that the rising edge is a hnear function of 
time in the vicinity of to- A similar technique has been 
utihzed with success by the VERITAS Cohaboration [S]. 
This assumption allows one to define the line by one sam- 
ple preceding and the other one following half of the peak 
sample Sp. These samples are labeled S- and 5*+ respec- 
tively. In this approach, the slope and intercept are given 

by 
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respectively, where the time i+ corresponds to the time of 
sample 5*+ and T is the sampling period, which is 4 ns for 
the two FADCs we studied. The characteristic time is now 
determined by 



tn = 



Sp/2 - bL 



a-L 



(7) 



In both approaches, we desire an algorithm that will 
provide a stable measurement of the characteristic pulse 
time independent of the pulse amplitude and other ran- 
dom variations in the pulse shape. In all cases the timing 
information will be used comparatively, that is, we are 
interested in comparing the characteristic time of some 
channel to another channel in the detector array. There- 
fore, constant, systematic offsets of to from the true time 
at which the signal reaches half maximum are irrelevant. 

4. Timing resolution measurements 

We tested the performance of the timing algorithms 
using an eight-bit single-channel and a twelve-bit multi- 
channel FADC. With the single-channel FADC we exam- 
ined digitization effects on the time resolution by effec- 
tively digitizing the same pulse twice with the FADC. Such 
an approach removes statistical variations in pulse shape 
and explores the true resolution of the algorithm. Using 
the multi-channel FADC we were able to recreate an envi- 
ronment comparable to that which we expect in the final 
detector. By examining channel-to-channel differences in 
the algorithm-determined characteristic time of a common 
light source we were able to measure final timing resolu- 
tions for individual channels and explore the dependence 
of these resolutions on pulse amplitude. 

4.I. Measurements with single- channel FADC 

A single-channel eight-bit prototype was used to mea- 
sure the timing resolution of the algorithms mentioned 
above. The light source in use was a small piece of scintil- 
lator, illuminated by a nitrogen laser LN300C. The light 
from the scintillator was propagated via optical fiber into 
a dark-box containing an FEU 84-3 and an XP-2020 pho- 
totube. The operating voltage for FEU 84-3 was set at 




100 150 200 250 300 

Sample Time [ns] 



350 



Figure 2: Typical FADC response for the primary and reflected pulse 
for two different delays. Comparing the reflected pulse between panel 
(a) and (b) demonstrates what is termed in the text a difference in 
"sampling phase." 



about 1600 V, and the fast XP-2020 PMT was used as a 
trigger. 

We devised a simple setup where the signal from a 
FEU 84-3 PMT was split, with one branch connected to 
the FADC and one end of the other branch left untermi- 
nated. The unterminated end produced a reflected pulse 
that was delayed, but still digitized within the same FADC 
buffer as the primary pulse. If one neglects distortion and 
attenuation introduced by the delay cable, this technique 
produces two identical pulses with fixed separation in time 
and removes contributions to the resolution from fluctua- 
tions in shape or transit time in the PMT. As long as 
the pulses are well separated in time, the algorithms de- 
scribed in Sec. [3] can be applied to obtain the difference 
in characteristic time of two pulses AtQ. Since the true 
delay is fixed, the standard deviation of Atg provides the 
quadrature-sum of the resolutions with which the charac- 
teristic times of the two pulses are determined. 

Figure [2] shows typical FADC responses for two differ- 
ent lengths of the open-end cable. Comparing the reflected 
pulse between panels (a) and (b) of Figure [2] demonstrates 
what we term a difference in "sampling phase." One can 
see that with Figure [2|b) the delay of the second pulse is 



such that the peak is hkely to fluctuate between two adja- 
cent samples, with the true peak being in the middle and 
slightly higher than the two. As both algorithms depend 
on obtaining a sample at the peak, one anticipates some 
added uncertainty when the primary and reflected pulses 
are sampled with different phases. By changing the length 
of the open end, one can control the delay between two 
pulses, modifying the probability that the sampling phase 
is different between the primary and reflected pulses, and 
explore the sensitivity of the timing resolution to the rel- 
ative alignment of the samples with the structure of the 
pulse. It is important to note that there is no synchro- 
nization between the FADC clock and the laser or trigger 
hardware, so it is impossible to control the absolute phase 
with which the pulses are sampled. 

The dependence of the time resolution on length of the 
delay cable can be seen in Fig. [3j where, for various delays 
of the reflected pulse, we plot the standard deviation of 
Aio, cr(Aio), obtained by applying the Gaussian (flUed cir- 
cles) and linear-interpolation (open circles) methods. The 
X-axis represents the mean of Aio modulo 4 ns. The best 
time resolution is obtained with a signal delay such that 
the peak region of both pulses are sampled with the same 
phase (as shown in Fig. I2ja)). Naively, one would expect 
this to occur only when the delay cable length is a multiple 
of 4 ns; however, bifurcated gaussian fits to both the pri- 
mary and reflected pulse indicate that the characteristic 
rise time, ai in Eq. [Tl of the reflected pulse is approx- 
imately 0.5 ns larger than the flrst due to dispersion in 
extra length of cable. This shifts the minimum of (T(Ato) 
from 4.0 ns (or 0.0 ns) to about 3.5 ns. The worst resolu- 
tion from both methods occurs when the delay is changed 
by approximately 2 ns with respect to this minimum, cor- 
responding to the case where both pulses are sampled with 
maximally different phases in the peak region (as shown in 
Fig.[2jb)). In this case, the resolution in time-difference is 
larger for the linear-interpolation method by about 70 ps. 

Assuming that CT(Aio) arises from equal contributions 
of the timing resolution of the primary and reflected peak, 
we can estimate the time resolution for the algorithm itself 
varies from 50 (40) to 110 (150) ps for the the Gaussian- 
transformation (linear-interpolation) methods. 

4-2. Measurements with multi-channel FADC 

Pulse-to-pulse and channel-to-channel variations in tran- 
sit time through the PMT and shape of the pulse will 
introduce fluctuations in the measured to - we call these 
fluctuations the "statistical" contribution to the resolution 
as they arise from random effects. We expect also an addi- 
tional "systematic" shift in the measured tg as a function 
of pulse height. This systematic shift must be corrected 
in order to compare to between two blocks in the array 
with different amounts of energy deposition. The studies 
presented in the previous section only explore the resolu- 
tion of the algorithm itself, hence a different approach is 
needed to examine these additional contributions to the 
timing resolution. 
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Figure 3: The value of (T(Ato) as a function of the mean delay (Atg) 
modulo 4 ns obtained by applying Gaussian (filled circles) and linear- 
interpolation (open circles) methods. The errors are statistical only. 



To better mimic the final detector configuration, we 
constructed an array of nine lead glass blocks optically cou- 
pled to PMTs which were powered by Cockcroft- Walton 
bases. All of the blocks in the array were simultaneously 
illuminated by a pulsed LED, which was mounted in front 
of a diffusion panel oriented about 30 cm from the face of 
the lead glass array. The LED was driven by a pulsing cir- 
cuit with a pulse rate of about 1 kHz. A capacitive circuit 
was used to shorten the duration of the applied electrical 
pulse to the LED to less than 5 ns. The entire setup was 
enclosed in a light tight box. One of the blocks in the 
experimental array was used as a trigger for the data ac- 
quisition system, while the remaining eight were connected 
to the FADC for data collection. 

Our goal was to extract the statistical and systematic 
contributions to the resolution as a function of pulse am- 
plitude. To reduce the pulse size, the incident light was 
attenuated. In order to have a set of control data with 
which to compare those blocks with varying pulse ampli- 
tudes, a screen was placed in front of only four blocks 
with the other four uncovered. The amplitude was varied 
by placing successively more screens in front of the ap- 
propriate blocks. The procedure was repeated using two 
different sets of control blocks. In order maintain a rela- 
tively constant transit time through the PMT throughout 
the study, the bias voltage on the PMTs remained con- 
stant at the values in Tab. [l] which roughly equalized the 
gains for all PMTs. (In the normal operating range of the 
FEU 84-3 PMT, variation in the bias voltage changes the 
transit time through the PMT by approximately 1 ns per 
100 V change in bias voltage.) 

As shown in Fig. |4J the difference in to for any two 
channels i and j for a set of events can be characterized 
by a Gaussian distribution with mean Aio,ij and stan- 



Table 1: PMT cathode voltage and pulse amplitude for various array 
elements. The typical amplitude of the pulses was about 500 mV 
when no light-attenuating screens were installed. 

Channel Voltage [kV] Amplitude [Channels] 



1 


-1.495 


1376 


2 


-1.545 


1398 


3 


-1.610 


1414 


4 


-1.585 


1405 


5 


-1.455 


1374 


6 


-1.500 


1401 


7 


-1.555 


1371 


8 


-1.515 


1439 



Table 2: Extracted tg standard deviations for the Gaussian trans- 
formation method (uq) and the linear interpolation method (tr^) at 
pulse amplitudes of about 500 mV. These measurements represent 
the dominant, statistical contribution to the resolution. 



Channel 


CTG [ns] 


(TL [ns] 


1 


0.28 ±0.01 


0.25 ±0.01 


2 


0.31 ±0.01 


0.28 ±0.01 


3 


0.29 ±0.01 


0.25 ±0.01 


4 


0.30 ±0.01 


0.28 ±0.01 


5 


0.27 ±0.01 


0.26 ±0.01 


6 


0.30 ±0.01 


0.28 ±0.01 


7 


0.27 ±0.01 


0.26 ±0.01 


8 


0.29 ±0.01 


0.28 ±0.01 
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Figure 4: This histogram of the difference in to in two channels shows 
that the distribution is essentially Gaussian in shape. Note that the 
two pulses do not arrive at exactly the same time. 



dard deviation aij . To obtain the mean value of the signal 
arrival time to^i and statistical contribution to the resolu- 
tion (Ti for a single channel, a simple set of equations was 
used. For n independent channels there are, in principle, 
n{n — l)/2 unique, non-trivial equations of the form 



Ato,^, 



= t( 



0,1 



o,i 



-I 



(8) 
(9) 



For a given configuration of the array, each pulse is digi- 
tized by the FADC and processed using the timing algo- 
rithms described in Sec. [3| The differences in to for suc- 
cessive channels are fit to a Gaussian shape, from which 
AtQij and aij are extracted. Collections of these fitted 
values are then utilized to study statistical fluctuations 
and systematic shifts in io as a function of pulse height, as 
discussed below. 

In order to obtain an accurate value of the statistical 
contribution to the resolution of a single channel as a func- 
tion of signal size, we first determine the resolutions for the 
four control channels by using a x^ fit. For these four chan- 
nels there are six independent values aij . It is possible to 



choose a set of four ai that minimizes a function given by 

4 4 / f^. _ Jcr? + a? \ 



X 




S{a.. 



(10) 



I]/ 



where S{aij) is the error in aij obtained from the Gaussian 
fit to the time- difference distribution. The minimization of 
this function yields values for the statistical contribution to 
the resolution. Repeating for the other set of four channels 
gives resolutions for all eight blocks in the array (at pulse 
amplitudes of about 500 mV) as shown in Tab. [2] 

Assuming the resolutions of the four control blocks are 
constant, the values of ai for the four blocks with atten- 
uated light input may be calculated using Eq. |9] The ex- 
periment is then repeated again with a different choice of 
control channels, allowing us to study all eight channels in 
the array. The process is repeated while varying quantity 
of light incident on the channels under study in order to 
map out the dependence of the statistical contribution to 
the resolutions as a function of pulse amplitude. The de- 
pendence is shown for the best and worst resolution chan- 
nels in Fig [5] for the Gaussian transformation and linear 
interpolation method. For relatively low pulse amplitudes 
the resolution is still better than 1 ns and appears to be 
in the range of 200-400 ps for larger pulses. 

The calculation of statistical contribution to the reso- 
lutions depends upon the choice of control channel used to 
determine aij. Any one channel may be compared to four 
different control channels. As a cross check on our proce- 
dure we varied the control channel used in our extraction 
of the statistical contributions to the resolution. This vari- 
ation produced changes in measured CTj of less than 40 ps 
over the entire range of pulse amplitudes indicating that 
our procedure is robust. 



The statistical contributions to the resolution were fit 
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Dnfj 



(11) 
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Figure 5: The statistical contribution to the timing resolution of 
channels 1 (open circles) and 8 (filled circles) are plotted as a function 
of pulse amplitude for the Gaussian transformation (a) and linear in- 
terpolation (b) timing algorithms. These two channels represent the 
highest and lowest resolution of all elements in the array - measure- 
ments for the other six channels are between these extremes. 



where a and b are free parameters (as shown in Figure pi). 
The average value for a over all eight channels was deter- 
mined to be 117 ± 34 (114 ± 46) channels-ns while that 
for b is 0.203 ± 0.079 (0.155 ± 0.077) ns for the Gaus- 
sian transformation (linear interpolation) algorithm. At 
an amplitude of 500 mV the average statistical contribu- 
tion to the resolution is 0.29 ± 0.09 (0.24 ± 0.08) ns, while 
at an amplitude of 100 mV the resolution is 0.62 ± 0.16 
(0.56±0.18) ns for the Gaussian transformation (linear in- 
terpolation) method. The uncertainties of the results are 
derived from the standard deviation of the resolutions for 
all eight channels studied. Both algorithms provide com- 
parable resolution. 

The values of the average pulse arrival time io,i calcu- 
lated from the same data sets are shown in filled squares 
in Fig. |6] These data show a significant systematic shift in 
io,i as a function of pulse amplitude. In order to account 
for this shift, which would otherwise cause instability in 
the measurement of to under variations of the pulse am- 
plitude, the values of to for each channel were corrected by 
the procedure 



to,i ~^ toA + OiSp 



(12) 



where a is a single parameter, apparently characteristic of 
the algorithm and FEU 84-3 pulse shape. In this study, a 
was calculated to be about —7 x 10^"* ns/ADC counts, in- 
dependent of channel and also digitization algorithm. This 
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Figure 6: The uncorrected to values (filled squares) as a function 
of signal amplitude calculated from the linear interpolation method 
for channels 1 (a) and 8 (b). The values after the correction given 
by Eq. |12| are shown as open squares. Other channels are consis- 
tent with these, and similar behavior is observed for the Gaussian 
transformation algorithm. 



correction, which is applied to the to the measured charac- 
teristic time, dramatically reduces the dependence of the 
characteristic time on pulse amplitude as shown in Fig. [6] 
and, consequently, reduces the systematic contribution to 
the timing resolution. 

The total contribution to the timing resolution from 
statistical and residual systematic fluctuations can be char- 
acterized as 



aiSp) 




(13) 



where a and b are determined as discussed above and c 
is the error due to the residual systematic contribution to 
the resolution. By computing the standard deviation of 
the corrected to for all channels over a range of pulse am- 
plitudes (as is shown for two channels in Fig. l6]) the value 
for c was determined to be 0.05 ns, a minor contribution 
when added in quadrature of the statistical contributions 
given by a and b. 



5. Conclusions and discussion 

We have presented two signal processing algorithms for 
converting flash ADC data into a measurement of the ar- 
rival time of an FEU 84-3 PMT pulse. One algorithm 
utilizes the fact that the leading edge of the pulse can 
be described by a Gaussian, while the other uses a linear 
approximation to the leading edge in the vicinity of the 
time at which the pulse has reached half of its maximum 



amplitude. We have measured the resolution of these algo- 
rithms and also the anticipated resolution with which we 
can measure a characteristic arrival time of a PMT pulse 
in the final GlueX detector configuration. 

By devising a technique to delay and digitize the same 
pulse twice with a single-channel eight-bit FADC we stud- 
ied the capabilities of the timing algorithms, independent 
of variations in pulse shape or transit time through the 
PMT. At a pulse height around 1.0 V, corresponding to 
roughly 85% of the full scale of an eight-bit FADC, we 
determine that timing resolution of the Gaussian transfor- 
mation algorithm is in the range of 50-110 ps, while the res- 
olution for the linear interpolation algorithm is 40-150 ps. 
The variation in resolutions is attributable to variations 
between the true time of the signal peak and the phase 
of the sampling clock. The time resolution that can be 
obtained by applying the two methods is only a fraction 
of the 4 ns time between samples. 

The relevant timing resolution for a single channel in 
the GlueX calorimeter includes contributions from the tim- 
ing algorithm, statistical contributions from pulse-to-pulse 
variations in shape and PMT transit time, and system- 
atic contributions arising from dependence of measured 
arrival time on pulse height. A nine-channel array was 
constructed to study these effects. The statistical fluc- 
tuations in arrival time, which depend on pulse ampli- 
tude, dominate the timing resolution and were measured 
as a function of pulse amplitude for both algorithms. The 
systematic shift in measured time as a function of pulse 
amplitude appears to be linear and consistent across all 
channels for both of the timing algorithms. This allows 
for a simple, single-parameter correction to be performed, 
which reduces the systematic contribution to the timing 
resolution to about 50 ps, making it nearly negligible. 
The linear interpolation algorithm performs slightly better 
than the Gaussian interpolation algorithm, although ei- 
ther would meet the resolution demands of the GlueX for- 
ward calorimeter. Assuming the use of the linear interpola- 
tion algorithm, the estimated final timing resolution antic- 
ipated in the detector for FEU 84-3 PMT pulses with am- 
plitudes of 100 mV (500 mV) is 0.57±0.18 (0.24±0.08) ns, 
where the central values include the quadrature-sum of 
statistical and systematic contributions to the resolution. 

Since the algorithms rely on sampling the pulse on its 
leading edge and near the peak, we anticipate that mea- 
surements of the characteristic pulse time will remain ro- 
bust in a high rate environment assuming PMT pulses 
do not overlap. The digital signal processing electronics 
could further be modified to apply the algorithm to mul- 
tiple detected peaks in a single flash ADC buffer. The full 
FEU 84-3 pulse flts inside of a 40 ns window. The 1-2% 
of FCAL modules that are closest to the beam axis have 
expected background rates of > 1 MHz. For the highest 
rate modules the probability of pulse overlap may reach 
10%, resulting in some small background contamination. 
For the remaining 98% of the modules, the probability of 
overlap is expected to be negligible. 



Both algorithms presented provide comparable timing 
resolution that is suitable for our application. The differ- 
ence in the two algorithms is in the implementation. The 
Gaussian transformation algorithm relies on more compli- 
cated mathematical functions; however, has the advantage 
that it always uses the peak sample and the two samples 
immediately preceding the peak. The linear interpolation 
algorithm, while mathematically more simple, requires a 
search and comparison of all samples on the leading edge 
to find the sample immediately before and after the point 
at which the signal level crosses half the maximum value. 
The exact choice of algorithm for the GlueX application 
will depend on which technique can be most efficiently im- 
plemented in flash ADC FPGA. 
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